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Regimes of Turbulent Combustion 



Ri 


Turbulent Reactive Flows edited by P.A Libby and F A Williams (1994) 


Regimes of Premixed Turbulent 
Combustion 



Turbulent Reactive Flows edited by P A. Libby and F A Williams H994) 




Regimes of Non-Premixed Turbulent 
Combustion 



Turbuient Reactive Flows edited by P A Libby and F A. Williams (1994) 


Chemical Closure Models 

(1) Laminar Chemistry 

< pWj >=pw i (Y i ,T) 


(2) Fast Chemistry 

i_~ a 2 v e (f) 

< pWi >=--p Zf ___ 


(3) Flamelet model 

< P w i >= JJ PWj (T|,X f )P f Xf (Tl,8 f )dnd£ f 


(4) Assumed PDF: 

<pw, > = J * * / P w i ( < l > i)’P^d(j) 1 d(j) 2 ..d(j) n 

Assumed the shape of P<|). 

(5) Scalar PDF method: 

Solve for P([) directly. 

(6) Conditional Moment Closure (CMC) 

<pw_ >=J<pw ; It) >■ P f (Tpdri 
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Flamelet Model: 69%H2+31%CH4 
Turbulent Jet Flame. Rev.=l 0.000 



VnuKK, « a! Nitnc Oxide Formation and Differential Diffusion m a Turbulent Methane- 
Hydrogen Diffusion FlamcT 24th Symposium(lmcmalionaJ) on Combustion/Tbc Combustion 


Flamelet Model: 69%H2+3 1%CH4 
Turbulent Jet Flame Rev.= 10.000 



Vranos. ct al Nitnc Oxide Formauon and Differential Diffusion in a Turbulent Methane- 
Hydrogen Diffusion Flame." 24th Symposiumf International) on Combusuon/Thc Combustion 
Institute. 1992/pp 377*384 
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extinction properties 


Advanced Flamelet Approach 


Local models 



Global model 


Conditional Moment Closure (CMC) 

Definition: 

<Y i lr|>^<Y (x,t)lf(x,t) = r) > 


Equation: 


, a < y, i-n > „ V1 v {<p U 'y'iii>p ! (Ti)} 

< plTl > 1 H < pulT| > V < Y- It) > H 


3t 


=< pwjTI > + < pDjVf • Vf iTl > 


P,(T1) 


a 2 < y : iti> 
an 2 


Modeling: 

< W| Iti >= Wj (< Tlr) >,< Yj lt| 

< pD^f • Vflr) >=< pD,Vf • Vf >= ^px, 

<pulT|>=pu 

< pu'y'lri >= 0 

< plTl >= p(< Y, lT| >,< Tit) >) 
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.Local strain rate, , 
'equivalence ratio... 


















Conditional Moment Closure (CMC) 

Applications: 

• Incorporated into existing moment 
closure CFD codes for complex 
geometry flows 

• Realistic Chemistry - Detailed or 
reduced 

Research issues: 

• Modeling of conditional statistics 

• Preferential diffusion 

• Parallel computing algorithms 


Probability Density Function (PDF) 

Applications: 

• NOx from methane jet flames with 
reduced chemistry 

• Sooting flames 

• 2-D flows 

Research Topics: 

• Mixing model 

• Extension to droplet spray & particle 
laden flows 

• Preferential diffusion 

• Efficient stochastic algorithm 

• Construction of chemical tables 

• Parallel computing - 3D Flows or 2D 
flows with complex chemistry 
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Mixing Models for PDF Methods 


- Modified Curl’s Model (stochastic) 
92 


k 

- I 


a=tp=l 3v a a V\ aP 


3 I*= v )/L ( mm 


1 1 » 


dy'd\i" 

mix xj/ty" 

<)> 0 <t> J 



- IEM Qnteraction-by-Exchange-with-the-Mean) Model 
(deterministic) 


k 

I 


Mixing Frequency 


10 


* 


' 2x . dxi!" 
mix Y a 




(0 = 
mix 


Modified Curls Model 



n 


IEM Model 



PaSR: H2/NOx Detailed Chemistry <j)=l x=lms 
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Comparison of Predicted and Measured 
H20 Mass Fractions 

Turbulent Nonpremixed Jet Flames 


Fuel: (CO/H2/N2: 0.30/0.10/0.6 ) 



Experimental Evidence of Preferential Diffusion 
in Turbulent Jet Flames 
(Fuel: 36%H 2 +64% C0 2 ) 




“Differential Molecular Diffusion in Reacting and Nonreacttnc Turbulent Jets of H2/C02 
mixing with Air.” L L Smith Ph D Thesis. University of California at Berkeley (1994) 


229 






Computation of Turbulent Reacting Flows 
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INTRODUCTION TO TURBULENCE SUBPROGRAM 
T.-H. Shih and J. Zhu 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 
Cleveland, Ohio 


OBJECTIVES 


• A means for CMOTT to interact with industry 

• A vehicle for technology transfer to industry 


CONCEPT OF TURBULENCE MODULE 
Exact CFD equations: 


d p u > . 9 r / dUi 9Uj 2 dU k dP 

Reynolds stresses will be recast ed as: 


= dUi . dUj 2 du k aUt ou 

pu t u j _ ^ r (^— — — -——Sij) + [-pu t Uj - ^z r (- — — — - 

dxj dx 1 


dxj dx t 3 dxk ' 


k 2 

VT = C»- 




CFD equations become: 


DpUi __ d 
Dt dx , 




k dxj dx x 3 dx k 


dxj dxj 


The task of turbulence module: Provide & and T 


CO I to 



• Turbulence Module: 

<> Input: Ui, p and p ... from the mean flow solver 
0 Output: 


fc 2 r Dk _ 

Mr - W — [— — ■ ■■, 


Dt 


De 

Dt 


= ...} 


m ,dUi , auj 2 du k , . 

T tJ - puiUj Hr( dx , + dXt zdxk * j ) 


<C> Models for p v~v- 

- One- and two-equation eddy viscosity models 

- Reynolds stress algebraic equation models 

- Reynolds stress transport equation models 
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Module with CMOTT research code (incompressible) 

• CFD equations in CMOTT research code: 


DpUi 

Dt 


d .. s .dUi dUj,, d „ 


dP_ 

dxi 


• Turbulence module: provide n T and T t] 

0 Built-in models without wall function: 

- Launder- Sharma and Chien k - e models 

- CMOTT k - e model 

0 Built-in models with wall function: 

- k - w model, standard k - e model 

- CMOTT k - e model 

- CMOTT Reynolds stress algebraic equation model 


Module with NPARC code 

• CFD equations in NPARC code: 


DpUi 

Dt 




Mr)( 


dU j dUj 

dxj dx t 


2 m s , 

3 dx k <l] 


dP 

dx 7 


• Turbulence module (present time): provide isotropic n T 
0 Build-in models without wall function: 

- Baldwin-Lomax model and Chien k-e model 

- CMOTT k-e model 
<> Further development: 

- Models with wall function 

- Reynolds stress algebraic equation models 

- Reynolds stress transport equation models 
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Joint Program with Industry 
on Turbulence Module 

• For those who want to use the available modules: 

0 Need interface program for particular industry codes 
- Grid informations, Boundary treatment, . . 

• For those who want a module for their own codes: 

0 Need modules exclusively for particular industry codes 

• Maintain and update the turbulence modules along with 
model development. 
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DESCRIPTION OF TURBULENCE SUB-PROGRAM 


J. Zhu 

Institute for Computational Mechanics in Propulsion 
NASA Lewis Research Center 
Cleveland, Ohio 


General Transport Equations 

- D i4>) = rJ^Sf 

• Non-dimensional form O, m fi/Re, m/Re ) 

• Conservative form 

• Cartesian velocity components 

1. Easy to transform (chain rule) 

2. No curvature terms 


235 



Discretization 


• Finite-volume method 

• Source term 

= S 1 + s 2 4 >, Si > o and S 2 < 0 

• Transient term 

1. lst-order fully implicit scheme 

2. 2nd-order three-level fully implicit scheme 

• Diffusion term 

Standard central differencing scheme 


• Convection term: HLPA scheme 
(Hybrid Linear/Parabolic Approximation) 

, / , L \i t _ - 4>ww 

4>W = <t>W + 7 (0C - (PwJVWi <PW — ( f )c _ ( p ww 


7 = 


1 

0 


if \4>w ~ 0 .5 1 < 0.5 
otherwise 


- Second-order accurate ■„ 

, J 1 I , 

- Bounded (non-oscillatory) ™ * 1 M e 1 

- Diagonally dominant matrix 

•SS 


EE 
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Solution Procedure 


• Non-delta form 

Positiveness (<f> > 0 but A0 may < 0) 

Simple linearization 

• Algebraic equations 

A c 4>c = A w 4>w + A e 4>e + As4>s + A n 4>n + S 
A's , S > 0 

• Decoupled solution 

• Alternating direction TDMA solver 


Boundary Conditions 

• Inflow: <p specified 

• Outflow: Fully-developed condition 

• Symmetry: d<fi/dn = 0 

• Wall: 

1. Low-Reynolds number turbulence models 

2. Standard wall-function approach 
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Sub-Programs 

• NPARC2D version 

Plane or axisymmetric, without swirling 

Compressible 

Non-vectorized 

• FAST2D version 

Plane or axisymmetric, with or without swirling 

Incompressible 

Vectorized 


NPARC2D Version 

9 ; I I ■ : i 1 : l 

[ 1 i r i ! j i i 

• Grid arrangement \ j { | I I | I j ‘-*--1 

5 : | • • ' ; ■ ' ! ; 1 

Control volume centers 1 3 j M ■ ] 1 : | f }-■ 

K , 1 M | ! T ■ - ' ! - 

Boundary nodes , * ? « 

J 

Embedded bodies 





zzzzzzzzzzzzz 
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J-Patches 
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• Input from the main code 

1 . Geometric quantities: x, y , €y, Px, Py, J 

2. Flow variables: p, J~ 1 p , J~ 1 pU , J -1 pV, J _1 E 

3. Patch control: 5x2 parameters 

4. Boundary conditions: 7x2 parameters 

• Output 

1. To the main code: pt 

2. For post-processing: K , e, ^ 


FAST2D Version 


• Grid arrangement 
CV centers 

B1 

Boundary nodes 
Embedded bodies 
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• Vectorization 
Single-index: 
ii=i+(j-l)ni 
0(iJ)=0(ii) 

0( i + 1 ,j ) =0(i i + 1 ) 
0(iJ-l)=0(ii-ni) 
Control parameter: 



KBLK = 


1 

0 


for computational nodes 
otherwise 


0=KBLK-0 c +(1-KBLK)0 6 


• Input from the main code 

1. Geometric quantities: x , y, x ^ x v , y^, y v , J 

2. Flow variables: p, p, U, V, W, C w , C s 

3. Vectorization parameters 

4. Boundary parameters 

• Output 

1. To the main code: pt, 

2. For post-processing: K, e, y+, y n , 
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OVERVIEW OF PROBABILITY DENSITY FUNCTION (PDF) MODELING AT LeRC 


D.R. Reddy 

Internal Fluid Mechanics Division 
NASA Lewis Research Center 
Cleveland, Ohio 


OBJECTIVE 

Accurately model the effect of turbulence on 
chemical reactions in a fluid flow 


APPROACH 

Use Probability Density Function (PDF) model - 
Express dependent variables as functions 
representing statistically realizable events 


POSSIBLE MODELING STRATEGIES 

1 . Evolution PDF - solve for function 

a. Joint PDF for velocities and chemical species 

b. Joint PDF for only chemical species 
& energy 

2. Assumed PDF - function prescribed 
Limited range of applicability - 

reaction time < < or > > turbulence time scale 
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CURRENT APPROACH 


• Develop evolution PDF model for compressible 
reacting flows & extend to spray combustion 

• Solve for joint PDF for species and energy 
using Monte-Carlo technique 

• Couple with conventional CFD codes 


AREAS OF IMPACT 

• NOx Prediction - HSCT and AST application 

• Spray combustion - swirling turb. reacting flows 

• Scramjet flow path analysis 

• Ignition kinetics - prediction of blow-off, etc. 

• Combustion instability studies 
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CODE FEATURES 


• Modular - can be coupled with any CFD code 

• Applicable for compressible flows with 
discontinuities 

• Monte-Carlo solver for generalized curvilinear 
coordinate system 

• Easily adaptable for parallel computation 
(currently under progress) 


CURRENT STATUS 

• 2-D and axisymmetric version released 

(default H2-air chemistry - 5 species) 

- parallel version to be released 

• 3-D version demonstrated for supersonic 
combustion (jet in cross flow) 

- validation planned for HSCT-type configurations 

• General chemistry (CHEMKIN) 

- Hydrocarbon spray combustion case currently 
under study 

• CFD codes used - RPLUS, ALLSPD, & SIMPLE-type 
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FUTURE PLANS 


• Further application/validation of 3-D model 

• Improved closure models - mixing and turbulence 
(use available DNS data) 

• Parallel processing - workstation clusters 

• Unsteady applications - long-term 

• Extend scope of impact 
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